Measurement of the short-range attractive force between Ge plates using a torsion 

balance 
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We have measured the short-range attractive force between crystalline Ge plates, and found 
contributions from both the Casimir force and an electrical force possibly generated by surface 
patch potentials. Using a model of surface patch effects that generates an additional force due to a 
distance dependence of the apparent contact potential, the electrical force was parameterized using 
data at distances where the Casimir force is relatively small. Extrapolating this model, to provide 
a correction to the measured force at distances less than 5 /im, shows a residual force that is in 
agreement, within experimental uncertainty, with five models that have been used to calculate the 
Casimir force. 

PACS numbers: 12.20.Fv, ll.lO.Wx, 73.40.Cg, 04.80.Cc 



Introduction. - The Casimir force has been a subject of 
great interest, both theoretically and experimentally, be- 
cause it is a macroscopic manifestation of quantum vac- 
uum effects P, 0, Q ' it '^^^ have significant effects in 
nanomechanical systems. Despite a number of success- 
ful measurements celebrated over the last decade [3, S] , 
early investigations of short-range forces 0, 01 report the 
possible systematic effects due to residual electrostatic 
forces. In particular, the observed variation in the effec- 
tive contact potential, recently reported in and later 
confirmed in Q , presents a problem of fundamental im- 
portance when setting limits to predicted submicron cor- 
rections to Newtonian gravity in a Casimir force measure- 
ment [loj . The optical response of a particular sample 



under study must also be carefully considered [ll|, as 
the accuracy of data on the optical properties of mate- 
rials typically limits calculational accuracy to no better 
than 5%. In principle, both electric and optical studies 
of a given sample are subject to a combination of various 
surface effects of electric origin, and it is important to 
understand these issues in order to accurately character- 
ize fundamental interactions, such as the Casimir force 
and non-Newtonian gravity. 

Our torsion balance set-up.- In this Letter, we present 
results of force measurements between crystalline Ge 
plates [12| in a sphere-plane geometry. Our apparatus, 
shown schematically in Fig. [l] is based on the design 



scribed in 0, 



presented in 13 1 and improves on the apparatus de- 



14l |. On one side of a torsion pendulum 



a flat Ge plate is mounted, and approached by a Gc 
plate with a spherical surface, with radius of curvature 
R = (15.10 ±0.05) cm, mounted on a Thorlabs T25 XYZ 
motion stage (8 nm resolution). When a force exists be- 
tween these plates, the torsion body rotates and thereby 
generates an imbalance in capacitance on the other side 
of the pendulum, which carries a flat plate, situated in 
between two fixed "compensator plates", that are at- 
tached to the support frame. An AC voltage is applied 



to the compensator plates, and the capacitance imbal- 
ance creates an AC voltage that is amplified and sent 
to a phase sensitive detector (PSD), providing an error 
signal to a proportional-integral-differential (PID) feed- 
back circuit. A small correction voltage (Spin) is applied 
to the compensator plates keeping the system in equilib- 
rium. The correction voltage is added to a large con- 
stant voltage Vo{~ 9V) to linearize the restoring force, 
F oc (S'piD + Vb)^ ~ + 2VbS'piD- This correction volt- 
age provides a measure of the force between the Casimir 
plates and is recorded during the measurement. 

The measured signal S'piD has contributions from sev- 
eral sources: 

Spmid, V,) = 5Dc(d ^ ^) + Sr{d) + Said, K), (1) 

where Sdc is the force-free component of the signal at 
large distances, Sr is the residual signal due to distance- 
dependent forces, such as Casimir-Lifshitz force, and Sp 
is the signal due to the electrostatic force in response 
to an applied external voltage T4. For the sphere-plane 
geometry, this latter signal can be written in the prox- 
imity force approximation (PFA), valid when d ^ i?, as 
Sa{d, Va) = 7reoi?(V'a — Vni)'^/(3d, where /3 is a calibration 
factor that converts S'piD in units of voltage to the ac- 
tual units of force. The electrostatic signal is minimized 
(Sa = 0) when = Vm, and the electrostatic minimizing 
potential is then defined to be the contact potential 
between the plates. 

A range of plate voltages T4 is applied, and at a given 
separation the response SpiD is fitted to a parabola 



SpiD{d,V^) = So + k{V,~V,af 



(2) 



The first two terms in Eq. ([T]) are absorbed in Sq and 
represent the minimized signal when T4 = Vm- Repeat- 
ing the parabola measurements shown in Fig. 2a, se- 
quentially moving from the farthest to closest plate sep- 
arations, enables us to inspect the d dependence of the 
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FIG. 1: Experimental setup of torsion balance (top-view). A 
pendulum body of length 15 cm hangs from a tungsten wire 
connected to a motorized rotation stage via the pivot point 
that is mounted on a support frame. The wire diameter is 25 
/im, with length 2.5 cm, shorter than the previous experiment 
(66 cm) d in order to minimize effects of tilt of the appara- 
tus. At the bottom of the pendulum body (not shown in the 
figure) is a NdFeB magnet to damp the swinging modes of 
the pendulum at a natural frequency of 3 Hz. The mechani- 
cal assembly is covered by a glass bell jar (vacuum 5 x 10^^ 
torr) and is supported on a vibration isolation slab that has 
its foundation separate from the laboratory building. 



fitting parameters k{d), Vni(d), and So{d). The proce- 
dure outlined here was first implemented as a calibration 
routine in [ist and more recently in Q in an effort to 
detect a distance dependence of 

As the gap between the plates is reduced, the parabola 
curvature k rapidly increases as shown in Fig. ^}p. These 
curvature values are fitted to k{d) = a/d, the expected 
dependence for the plane-sphere geometry, where the ab- 
solute distance d = do — d^ is defined in terms of the 
asymptotic limit do and the relative distance dr recorded 
during a parabola measurement. The conversion factor f3 
is then obtained through a = 'ne^R/ f3. Obviously, a can 
be also used to determine the absolute distance through 
d = a/k, implying a significant correlation of a with do- 
Consistency between these two methods of distance de- 
termination reflects validity of the use of the \/d power 
law as implied by a value of Xq close to unity for our data 
set. Fig. [2t shows the electric potentials Vm at minima of 
the parabola curvatures plotted versus d, indicating the 
distance-dependent minimizing potential Vm{d)i a behav- 
ior that has been observed in other experiments 0, 0] . 

To see the trend in Vm(d) more clearly and to de- 
termine short-range forces with higher statistical accu- 
racy, we have repeated 200 times the experimental se- 
quence described in Fig. [21 yielding a total of 5800 data 
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FIG. 2: Description of procedure for force analysis, (a) A 
single parabola measurement at a given distance is acquired 
by sweeping a range of voltage differentials applied to the 
plates. The procedure is repeated at decremental distances 
from 150 /xm down to 500 nm, completing a single experi- 
mental run. The parabolic curves shown here represent only 
three distances, 100 ^m, 20 /im, and 1 yLxn for clarity, (b) Cur- 
vature coefficients of the parabola k versus relative distances 
were fitted to the 1/d electric force to provide the voltage- 
to-force conversion factor (3 as well as the absolute distance 
obtained from an asymptotic limit, (c) The force-minimizing 
potential Vm changes with distance, varying approximately 
as log d. (d) The residual force at the minimizing potential is 
plotted against distance, after subtracting the DC offset 5dc 
and multiplying it by (5. The maximum force gradient for 
feedback system stability is 5 uN/pm, limiting the minimum 
distance to 500 nm. 



points. Each group of five data points taken at a given 
fixed distance with varying applied potential are used 
to determine the three parabola parameters discussed 
above, in addition to the force and distance. The mean 
value of the calibration factor after analyzing all data is 
13 = (1.35 ±0.04) X 10-^ N/V. Both the asymptotic limit 
do, shown in FigOa, and the DC offset of the PID signals 
5dc drift slightly during a run. The uncertainty in po- 
sition is roughly 10 % at a given distance and about 50 
nm at the typical closest gap separation, consistent with 
the actuator minimum displacement of 40 nm. The DC 
offset drift has been corrected by monitoring 5piD be- 
fore and after each consecutive run and applying a linear 
correction. 

Varying minimizing potential.- An outstanding feature 
of our data is the distance variation of the applied volt- 
age Vin that minimizes the force, as clearly shown in Fig. 
[3l It must be recognized that this variation can lead 
to an extra force of electrical origin, as demonstrated in 
However, the model used in [3l assumes that the 
variation in the minimizing potential is due to a vary- 
ing contact potential, specifically modeled as a voltage 
source in series with the plates. The varying minimizing 



3 



30 



> 



I 27 

a 
S 
o 
& 

00 

■ S 24 
N 



21 





3,11 



4 r 



10 



10"^ 

Distance (m) 



10 



FIG. 3: The value of the force minimizing potential as a func- 
tion of plate separation. The red points are the average of the 
data. Our measurement reveals a slow rise of the minimizing 
potential as the plates approach each other, of order of 6 mV 
over 100 ^m. This variation in our data set shows a sim- 
ilar trend observed in a recent measurement @] where the 
variation of 6 mV over 1 /im is reported between An coated 
samples, significantly larger than the value found in our study. 



potential observed in our data is more likely due to large- 
scale gradients in the contact potential across the surface 
of the plates, due to. for example, polishing stresses or 
the curvature of the spherical surface plate changing the 
crystal plane orientation at the surface. Such variations 
have been observed for many materials [itI . , with typ- 
ical large scale fluctuations on the order of a few mV. The 
variation in the apparent contact potential is due to the 
effective averaging area changing as the curved and flat 
surfaces of the plates are brought together. A numerical 
analysis [isj of a wide range of surface potential varia- 
tions shows that the variation of Vni{d) leads to an elec- 
trostatic force of the form = 7ri?eo[Ki(rf) + 
where Vi is a constant offset parameter at large distances 
that can be determined from the experimental data. The 
origin of this effect is due to the plate curvature, together 
with large scale variations in the surface contact poten- 
tial. Note that although the parabola measurement min- 
imizes the electrostatic force across the plates, it does 
not necessarily nullify all the electric forces that possibly 
exist. 

Random small-scale patches. - In addition to large scale 
gradients in the surface potential, there can be small- 
scale (i.e., much smaller than the plate diameter) ran- 
dom fluctuations in the surface potential associated with 
strains, irregularities, and impurities. It is straightfor- 
ward to show that the electrostatic energy per unit area 
between two flat plates with random patch voltages is 



2 

rms 



g-2M 

duS{u/d)— — 2 J 

sinh (u) 



(3) 



where K-ms is the rms value of the random patch 
voltages. For simplicity, we have assumed isotropic 
patches with surface correlation functions (Vk^iVk'j) = 
V^^^^S{k)SijS{k — k'), where i,j = 1,2 denote the plates, 
and >S'(fc) is the unity-normalized spectral density. The 
residual electrostatic force between the sphere and the 
plane due to these patches can be obtained from PFA 
as Ff^{d) = 27ri?i?patch(rf)- For example, for random- 
voltage patches of radius A uniformly distributed on the 
surfaces, the spectral density is S{k) ss sinfcA/TT^/c^. 
It is easy to see that in the limit A 3> d the resid- 
ual patch force in the sphere-plane geometry scales as 

Ff^[d) = TrReoVl,JdM- 

Electrostatic residual force.- We fit the data of the 

residual force at the minimizing potential (Fig 2.d) with 
a force of electric origin Ff = F^ + Ff^ + Ff^ = F^ + 
7ri?eo{[Ki(rf) + ViY + K-msl/c'i where Fq is an offset pa- 
rameter at large distances. A least-squares fit of the ob- 
served force using data for d > 5 ^m (a regime where the 
Casimir force should be vanishingly small) and the mea- 
sured Vni{d), while leaving Fq, Vi, and K-ms as adjustable 
fit parameters, yields an excellent description of the ob- 
served large distance force 2l| , as shown plotted with the 
data in Fig. 4. Including data at shorter distances (d < 5 
fim) causes a significant fit deviation, indicating an in- 
terference with the actual Casimir force which is highly 
nonlinear at short distances. A similar long-range force 
has been previously observed in the measurement of van 
der Waals interaction and the corresponding correction is 
applied to the data based on work function anisotropics 
and their related patch charges Q. 

Casimir residual force.- We have tested our data for 
the presence of a residual Casimir force F^'^{d) be- 
tween the Ge plates, which we have computed in PFA 
from the plane-plane Casimir-Lifshitz energy, F^'^^{d) = 
2TTRE^p^{d). We have calculated the corresponding re- 
flection coefficients using five different theoretical models 
for the Ge plates 23|, l2J]: ideal dielectric; ideal dielec- 
tric + Drude conductivity corrections; ideal dielectric 
+ plasma conductivity corrections; quasi-static Debye- 
Hiickcl screening model; and charge drift model. Fig. 
4 (bottom) shows the experimental data for the resid- 
ual force after subtraction of the joint (contact potential 
and patch potential) electrostatic forces, and the theory 
curves for the Casimir-Lifshitz force between Ge plates 
at r = 300 K for these five theoretical models. The error 
bars take into account all statistical uncertainties (2-3 
%) as well as fitting uncertainties from the electrostatic 
force analysis (10%). Within experimental uncertainty, 
our data agrees well with the theoretical predictions for 
the Casimir force, but it is not of sufficient accuracy to 
distinguish among the different models for the Ge plates. 
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FIG. 4: (top) Experimental data from parabola measure- 
ments along with a solid line representing an average. The 
joint (contact potential and patch potential) electrostatic 
residuals is modeled by F^^ = Fq + 7r7?eo{[V;^(d) + Vif + 
K-ms}/'^! which is fit to our data points, yielding Fq — 



(-if ± 2) X iO"'" N, Vi = (-34 ± 3) mV, 



(6 ±2) 



mV, and Xo = 1-5. (bottom) Experimental data for the re- 
maining force after subtraction of electrostatic correction 
together with theoretical Casimir forces computed using five 
models described in |2^,[23l. for d < 5 fim is approximately 
unity with the corrected data for all of the models. 



High precision measurements of the Casimir force require 
a careful evaluation of the electrostatic effects considered 
in detail for the first time in our present study. 

Conclusions.- We have performed measurements of the 
short-range force between Ge plates in the sphere-plane 
geometry, and have observed that the potential Vm that 
minimizes the electrostatic force depends on the gap be- 
tween the plates. We have considered two contributions 
of electric origin present in the residual data for the 
force. The first contribution is due to large-scale vari- 
ations in the contact potential along the surface of the 
plates, that leads to the gap-dependent minimizing po- 
tential and, as a result, to an electrostatic force propor- 
tional to {Vin{d) + Vif /d. The second contribution can 
be modeled as arising from potential patches on the sur- 
faces that, in the case when they have typical sizes much 
smaller that the plate diameters and much larger than 
the plate separation, leads to a further electrostatic force 
proportional to V^^^/d. We have fitted our experimental 
data at large distances (d > 5 ^m, where the Casimir 
force is expected to be negligible) with these two elec- 
trostatic force effects, and found we could establish good 
agreement between our model and the experimental data. 
Furthermore, we have subtracted these forces from the 
data at short separations (d < 5 /im) and found a resid- 
ual force that is in agreement with the theoretical pre- 
dictions for the Casimir-Lifshitz force between Ge plates. 



Our measurements do not have enough accuracy to dis- 
tinguish between the different theoretical models used to 
characterize the optical properties of the Ge plates. Fu- 
ture measurements are deemed necessary in light of our 
discussion, in particular to better understand the phys- 
ical origins of the observed electrostatic forces. We are 
currently exploring similar surface effects in a pair of Au 
samples. We acknowledge support from Yale University 
for the construction of the experimental apparatus and 
data acquisition, and from Los Alamos LDRD program. 
We thank G. Klimchitskaya for useful discussions. 
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